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Vacuum polarization mediated by quark loops is susceptible to external electro- 
magnetic fields as well as to the QCD vacuum structure. Employing the stochastic 
vacuum model, we calculate the modification of the one-loop Euler-Heisenberg 
effective action due to stochastic color fields with the Fock-Schwinger technique. 
Our results indicate nonperturbative light quark contributions of the same order 
of magnitude as the usual QED terms. Various theoretical and experimental im- 
plications are discussed in this progress report. 



1 Introduction 

The physical vacuum has become a most important research topic with the 
advent of quantum field theories of the fundamental interactions, of the gauge 
field theory of the Standard Model in particular B Basic properties of matter, 
such as masses of elementary particles and ultimately the observed particle 
spectrum, are induced by the vacuum state and its (broken) symmetries, which 
in turn are determined by the underlying interactions. Therefore, it is of 
considerable interest to study in any conceivable way the dynamical vacuum 
and, especially, its altogether still unknown quantum wave functional. 

Our aim here is to perform a first step in studying the interference between 
the very different color-confining QCD and charge-screening QED vacua. Ma- 
jor efforts have been launched to investigate the 'melting' of the QCD vacuum 
and formation of a quark-gluon plasma in high-energy heavy-ion collisions, see 
Ref.EI and earlier references therein. However, relatively little attention has 
been paid so far to the mutual influence of QCD and QED on the vacuum 
state and its virtual excitations. Naturally, since quarks carry color as well as 
electric charges, both contribute to vacuum polarization. 

Following earlier theoretical studies of the vacuum induced nonlinear photon- 
photon interactional there are ongoing searches for the QED vacuum birefrin- 
gence effect in high-precision crossed laser/magnetic-field experiments.! Most 
recently, we proposed to evaluate the influence of the_QCD vacuum structure 
on the Euler-Heisenberg (EH) effective action of QED .EM This well-known one- 
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loop effective action has been the central quantity for the evaluation of vacuum 
polarization effects in external electromagnetic fields, mosLnotably Schwinger's 
calculation of pair production in a constant electric fieldE 

Our intention is to calculate the nontrivial correction to the EH effective 
action arising from QCD interactions and the QCD vacuum, in particular. In 
lack of a solution of the QCD vacuum problem from first principles, we resort 
to the stochastic vacuum model (SVM)Ja which successfully describes color 
confinement, the quarkonia phenomenology, and produces interesting results 
in high-energy scattering, see RefJj and references therein. This model allows 
us to calculate the one-loop QCD vacuum-to- vacuum amplitude in the presence 
of external electromagnetic fields. 

By means of the SVM, the action-weighted functional integration over 
gluon (gauge potential) fields, S° a = 1 , . . . , N% — 1 for SU(N C ), is replaced 
by an ensemble average over a Gaussian distribution of the color fields, G®, in- 
stead. The latter is characterized by a single nonvanishing vacuum correlation 
function, which has been computed in nonperturbative lattice gauge theory E3 
and is related to the vacuum correlator appearing in QCD sum rules. LU 

We apply the Fock-Schwinger technique in order to calculate the one-loop 
effective action of quarks interacting nonperturbatively with an ensemble of 
time independent homogeneous color fields as well as with an approximately 
constant electromagnetic background field. 

Similarly—as the colormagnetic instability of the Mantinyan-Savvidy vac- 
uum ansatz,E3 the SVM points to the importance of inhomogeneous and pre- 
sumably time dependent color fields, which produce a characteristic correlation 
length for the QCD vacuum. By the uncertainty principle, a typical vacuum 
polarization loop has an extent of I « At • c ss hc/(2mc 2 ), which amounts to 
l e w 200 fm for electrons and l u « 20 fm for the lightest up-quarks. Of course, 
the physically relevant value of l u must be expected to differ considerably from 
this estimate guided by perturbation theory. Results from lattice calculations 
indicate that important nonperturbative corrections have to be taken into ac- 
count. One may consider an estimate of a lower bound to be determined by a 
typical constituent quark mass of toq ~ 300 MeV instead, i.e. Iq rs 0.3 fm. 

Whereas all macroscopic (laboratory) electromagnetic fields are constant 
on the scale of l e , the QCD vacuum correlation length obtained from lattice 
calculations, A « 0.2. . . 0.4 fm, is of the same order of magnitude as Iq and, 
therefore, is expected to remain small on the scale of the light quark loops 113 

Therefore, truly space-time dependent color fields also have to be incorpo- 
rated in the future. However, their study needs a technically different approach 
still to be developed. 



2 



2 One-loop QCD/QED Effective Action 



We briefly rederive the one-loop effective action for fermions in color plus 
electromagnetic background fields, aiming at the case of a constant electro- 
magnetic field strength with stochastic gluon field fluctuations, in particu- 
lar. - The calculations are performed for Minkowski space with the metric 

— (1, — 1, — 1, — 1), any Lorentz or color indices occuring twice are to be 
summed over, and we choose units such that Ti = c = 1. 

We are interested in the vacuum-to-vacuum amplitude which determines 
the QCD effective action, T^, in the presence of an electromagnetic background 
potential, A^: 

exp{iT A ) = J V^VipVB exp U J d 4 x (L A , B - X -F 2 - ^G 2 )^ , (1) 

where L denotes the Dirac Lagrangian for the quarks coupled to the electro- 
magnetic and color gauge fields, A and B, respectively, L a .b = V>(«7 ■ D — m)ijj, 
with D = d — ieA — igB denoting the covariant derivative in the funda- 
mental representation of the color SU(N C ) group; eventually, one has to sum 
over the contributions of all fermions, i.e. with different charges and masses. 
The gauge field actions are written in terms of the field strength tensors, 
F^u = d^Av - duA^ and GJ„ = - d„B° + gf ahc B\Bl, respectively. 

We did not add gauge fixing terms, since the gluon gauge potentials will be 
treated in terms of stochastic background field strengths in this work. We re- 
mark that the vacuum contribution to Y a is selected by applying the Feynman 
boundary condition, which will be made explicit shortly. 

In order to proceed, we define the Gaussian stochastic ensemble and its 
correlators,!! which will be employed in the following evaluation of the effective 
action. We assume for simplicity that only approximately constant colormag- 
netic fields contribute which, however, fluctuate in direction and amplitude. 

As previously pointed out in Refs.tij in the context of QCD background 
field calculations, the commutator of two covariantly constant color field strength 
tensors vanishes, [G^,G pa ] = 0. Consequently, they can be parametrized as 
being proportional to the generators of the abelian Cartan (sub-)algebra of the 
SU(N C ) color group. 

We thus take into account SU(3) colormagnetic fields of the form B = 
i(B 3 A3 + B S A$), involving the diagonal Cell-Mann matrices A3 ; s. Then, their 
normalized probability distribution is defined by: 

P(gB a ) d 3 (gB a ) = ( ^ ) 3 / 2 exp( ^ gB ^ ) d 3 (gB a ) , (2) 

M9 2 B 2 )g 2(9 2 B 2 )g 
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for a — 3, 8 , with the width determined by the relevant correlator. 

In accordance with the SVM modelja we replace the functional integral 
J VB in Eq. ([!]) by the Gaussian ensemble average with the distribution of 
Eq. (|^), which we denote by (. . .)g from now on. 

In the absence of a calculation from first principles, the correlator in 
Minkowski space is obtained by analytical continuation of the Fourier trans= 
form of the euclidean one, which is calculated numerically in lattice studies.li3 
This procedure is discussed in detail in Ref.El. Since we consider only stochas- 
tic colormagnetic fields, we identify the widths of their distributions with the 
value of the SVM field correlator: 

(y^GDc = (9 2 B 3 ■ B 3 ) G + (g 2 B s ■ B 8 ) G = 2(g 2 B* ■ B 3 ) G , (3) 

where g denotes the renormalized QCD coupling. The last equality (isotropy 
in color space) allows us to relate the respective correlators to the physical 
gluon condensates 

(-0^)0 = (- : G a G% u :) = 0.024 ± 0.011 GeV 4 , (4) 

where the (running) strong coupling constant is a s ,y .g 2 ./4ir, and we cite the 
empirical value of the condensate (» (394MeV) 4 )Jl3'BllJ 

For a selfconsistent determination of the stochastic field ensemble, the vac- 
uum condensate in particular, the gluonic contribution to the effective action 
has to be calculated incorporating gluon fluctuations in the stochastic back- 
ground together with an appropriate gauge fixingJl3 Here we limit our attention 
to the fermionic contribution, which generalizes the familiar Euler-Heisenberg 
effective action of QED. 

Obviously, the usual EH effective action EIlI for leptons, Tgjj, is obtained 
from Eq. (|l]) by dropping completely the integration over the gluon field B and 
setting the corresponding terms of the classical action to zero. In order to 
evaluate Ta instead, employing the SVM model with the stochastic ensemble 
average defined in Eqs. (||)-(|]), we begin by integrating out the fermions: 

iT A =\n((&et{[ 1 -D + im][ 1 -d + im]- 1 }) G ) , (5) 

where we subtracted the contribution of the nonintcracting fermionic vacuum 
fluctuations for later convenience. From now on, the Feynman boundary con- 
dition is implemented by replacing m — > m — ie, which we often suppress. 
The above result, Eq. (^), is similar to the usual QCD effective action (for 
= 0), however, with the vacuum correlator appearing implicitly instead of 



4 



the gluon propagator. Next, interchanging the order of taking the logarithm 
and the ensemble average and using det{M} = expTrln{M}, the result is: 

*(T A ) 1 _ loop = (TT\n{l 1 -D + i m }[ 1 -d + i m }- 1 }) G , (6) 

with a trace over space, spin, and color. We remark that interchanging the 
logarithm with the ensemble average results in the 1-loop approximation here 
(suppressing this subscript henceforth). Since the dynamics of the Dirac field 
does not depend on the sign of the mass term, taking the average of both 
possibilities, the 1-loop action can be expressed conveniently, 

zT A = \ (Tr ln{[(7 ' D) 2 + m 2 ][( 7 • df + m 2 ]" 1 }) Q 

= \ (Trln{[II 2 -m 2 + \a, v {eF^ +gG^)}[P 2 - m 2 ]" 1 }^ , (7) 

where we introduced the kinetic and canonical momentum operators, IT M = 
iD^ = + eA^ + gB^ and P M = id^, respectively; the second equality fol- 
lows with the help of {7 M , 7„} = Ig^v and [7^,7^] = —2ia^ v , i.e. the (anti-) 
commutation relations of 7-matrices. 

Finally, using the integral representation of the logarithm which presents 
the starting point of the Fock-Schwinger proper time method, we obtain the 
familiar looking result :0 

/*°° A 1 

T- A = — (Tr{exp(i S [n 2 - m 2 + -a^eF^ + gG^)]) 

Jq L 

-exp(z S [P 2 -m 2 ])}) G , (8) 
with the mass (— » m — ie) incorporating the Feynman boundary condition. 



3 Fock-Schwinger Technique 

For arbitrarily varying background fields it is unknown how to evaluate Eq. (|J), 
even for the QED case. Presently, we work with the simplified ensemble of color 
background fields introduced in the previous section, which have a fluctuating 
amplitude with respect to space and (color) space direction but are covariantly 
constant otherwise. We also incorporate the approximately constant electro- 
magnetic field to all orders, expanding the results only in the end. 

We proceed by relating the first exponential in Eq. (|J) to an unitary evo- 
lution operator, U(s), depending on the proper time variable s: 

U(s) = exp(-isH) , H = Ii 2 -m 2 + ^a fil/ (eF»» + gG» v ) , (9) 
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where H plays the role of a fictitious Hamiltonian. Following the Fock-Schwinger 
strategy as described in Refflt£l, the aim is to obtain the coordinate space ma- 
trix elements of H. Using these, the equation of motion for the evolution 
operator can be converted into an ordinary differential equation. Its solution 
provides the coordinate space matrix elements of U(s), which are sufficient to 
calculate the trace and stochastic average in Eq. (||). 

To begin with, introducing the Heisenberg operators n(s) = U^(s)WJ(s) 
and x(s) — [/^ (s)x£/(s), we obtain the related Ehrcnfcst equations of motion: 

d sXfl = i[H,x^] =-2n M , (10) 

s n M =i[#,n M ] = 2{eF^+gG^)Ii x +igJ^ + ^g[D^a-G] (11) 
w 2(eF fiX + 9 G M )n A , (12) 

where the color current, = [D , G^x], is introduced to indicate the phys- 
ical meaning of this term; similarly, the last term contributes the spin-color 
coupling here, with a ■ G = a^G^ , while the first term is related to the elec- 
tromagnetic and color Lorentz forces. In order to arrive at Eqs. ([To]) and ([TT1), 
we made use of the coordinate representation of the operators II, x introduced 
after Eq. ([?]) and treated as constant. The last equality, Eq. (|lj), presents 
the approximation for covariantly constant fields, i.e. our starting point here. 
It is convenient to introduce a combined field strength tensor, 

= eF^+gG a ^ , (13) 

which is a matrix in Lorentz and color indices, with t a , a = 1, . . . N% — 1 
denoting the generators of SU(N C ) in the fundamental representation. Then, 
the solutions of Eqs. ([l(]) and ([12]) are easily obtained: 

n(s) = ex P (2^ s )n(o) = — r^jTT—i M*) ~ x ^ - ( 14 ) 

exp(— zj-sj — 1 

^) = 1 ~ ex P (2 - Fs) n(o) + ^(o) , (is) 

where we employed Eq. ( |l5|) to eliminate 11(0) from Eq. (|l4|). 

With the help of Eq. (|15[) and the basic commutator [11^(0), 2^(0)] = ig^ 
one derives: 

[x,{r),x v {r')] = -2i (cxp^r - /]) = T ' D ) . (16) 
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Using this, the Hamiltonian of Eqs. (|^) can be written in-.time-ordered form 
and, then, its coordinate space matrix elements evaluated ]ZrlL£l We obtain: 



H(x',x; s) — ~[x' — x] 5 i x ' ~ x ] ~ o^r_L [.Fcoth(.Fs)] — m 2 + —a ■ T 

4 sinh (Ts) 2 



lr . , T 2 , , , i r _ . „ , 1 

7 ■ 

(17) 

where the trace refers to the Lorentz indices. 

Next, we turn to the equation of motion for the proper time evolution oper- 
ator U(s), which follows from Eq. ([)]). Interchanging between the Schrodinger 
and Heisenberg picture, we obtain in the coordinate representation: 

d s U(x',x;s) = d s (x'\U(s)\x) = a,(ar'(a)|aj(0)) = -i(x'\H U(s)\x) 

= -i{x'(s)\H(s)\x(Q)) = -iH(x'(s),x(0);s){x'(s)\x(0)) 

= -iH(x',x;s)U(x',x;s) , (18) 

with H{x' ', x; s) from Eq. (^). This presents an ordinary first order differential 
equation for U(x ,x;s) as a function of s, which can be integrated directly. 
Since all matrices involved are considered as constants at this point, no ordering 
prescription for the resulting exponential is needed: 



U(x',x;s) = C[x',x)exp(-i J dr H {x' , x; r)) . (19) 
The function C(x',x) incorporates the usual QED Schwingcr string EE 

C(x',x) = Cexp(ie / dz^A^z)) , z(£) = x + £(y - x) , 0<£<1 , (20) 

J X 

as well as a normalization constant C. Furthermore, we calculate: 

i , , 1 sinh(jFs) 



f s i If 

i / dr H(x',x;t) = -[x' — x]T coih(T s)[x' — x] — -tr L [ln(- 



+ ln(s~ 2 ) + im 2 s - ■ Ts , (21) 

where we separated the second logarithm for later convenience. 

We remark that the Hamiltonian, Eq. (IHf), is covariant under the global 
SU(N C ) gauge transformations admitted here and invariant under arbitrary 
electromagnetic gauge transformations. Therefore, the evolution operator, 
Eq. (|T9|), requires the additional electromagnetic string factor, in order to be 
properly covariant under either gauge transformations. 



Finally, the normalization constant C takes the boundary condition (or- 
thogonality and normalization of coordinate eigenstates) into account, 



\imU(x',x; s) = lim(x'(s)|x(0)) = 6 4 (x' - x) ■ l sc , (22) 

s— *0 s— >0 

with a unit matrix for spin and color on the right-hand side. The normalization 
constant is calculated using Eqs. ([l7|)-(|2^). The result is: 

Furthermore, Feynman's m — > m — it provides the convergence factor in 
Eq. jl9|), thus implementing the asymptotic condition U(x' , x;s—> — oo) = 0. 

This completes our derivation of the (matrix elements of the) evolution 
operator U. We note in passing that these results immediately yield the prop- 
agatoc,p. the considered combination of color and electromagnetic background 
fields,0H3 which may be useful for other purposes. 



4 QCD-modified Euler-Heisenberg Lagrangian 

Here we employ our assumptions about the nature of the stochastic color fields, 
in order to calculate the effective action, Eq. (|8|), using the results of the pre- 
vious section, Eqs. (|l9|)-(p3|) in particular. 

We begin by evaluating the traces over spin and space in: 

tr S , x U\s) = ^— ^ J d 4 xtr s {cxp(--tr L [ln( ^- > )]-im 2 s + -<j ■ Ts)} . 

(24) 

We recall that the covariantly constant color fields are parametrized propor- 
tional to the generators of the abelian subgroup of SU(N C ), as discussed in the 
context of Eqs. (|))~(U), i.e. T is a diagonal color matrix here. Then, the trace 
over spin is calculated similarly as in the QED easel] 

tr s exp(^cr ■ Ts) = 2(cos[^v/Xi + il 2 ] + cos[|-y/2i ~ ^2]) > ( 25 ) 
where two Lorentz invariants of the fields enter, 

J 1 = 2^ tlv ^ v , l 2 = e al3lS T a0 ^ 5 . (26) 
The remaining exponential factor in Eq. ( p4| ) can also be simplified: 

1 ,sinh(J r s) —is 2 ,s , s , 

exp --tr L [ln( )]) = — — T 2 (cos [-y/li + i2" 2 - cos ~ ^2 , 

2 J- s 4 2 2 

(27) 
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again quite similarly as in the QED case. El 

Combining Eqs. (p4))-(p7|) and using these in Eq. (|J), we obtain the un- 
renormalized 1-loop effective Lagrangian (T A = J d 4 x C A ) : 

C- A = (-l/8^ 2 ) • (28) 



where the color trace and the stochastic average are still left to be done. IfUwx 
omit these and set all color fields to zero, the usual QED result is recovered.LHI^ 



4-1 Renormalization 

Since QCD and QED are both renormalizable, we are guaranteed that the 
ultraviolet (s — > 0) divergences contained in C A can be absorbed by renor- 
malization, of the charges and fields in particular. At present we are mostly 
interested to demonstrate how QCD modifies the usual Euler-Heisenberg ef- 
fective Lagrangian. Therefore, it is convenient to subtract from C A the pure 
QCD contribution, Cq . 



It is calculated by setting e = in Eq. (28). Then, for colormagnetic 



vacuum fields, we obtain the renormalized result: 

C o = -g^< ^c{g 2 B ■ B / ^e-* a [<rooth(aO - 1 - -x 2 }} ) G , (29) 



where z = m 2 / y g 2 B ■ B is a diagonal color matrix. We remark that here we 
rotated the integration s — > — ix in the complex plane, as compared to Eq. (psh; 
it is obvious now that the integral increases strongly with decreasing m. The 
corresponding QED integral has been calculated analytically.tll Using this, we 
obtain in the limit of strong fields (z <C 1): 



C ° = 24^ ( tYc{g2§ ' B H™ 2 /V.9 2 5 ■ B) + 0(z°)]} ) G , (30) 

which is the appropriate limit for the light quarks, given the large value of the 
gluon condensate, cf. Eq. (|J). 

After the QCD renormalization, we have to replace gG^ = gnG^u also 
in the remaining terms, i.e. by the renormalized quantities; we drop the sub- 
script R henceforth. The resulting QCD- subtracted effective Lagrangian is: 
C A = C A - Co , cf. Eqs. ©-© . 
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In order to proceed, we introduce some useful abbreviations: 



a=^{2F^F» v -ie^gF^F^) , d = a* , (31) 
b = ^(2F^G^ - ie a0lS F a ^ s ) , ~e = b* , (32) 

~c=^(2G^G» v -ie aPlS G al3 G~< s ) , f = c* , (33) 

where i 7 ^^ and G Mt , are understood as a unit and a diagonal color matrix, 
respectively. Then, the field dependent factor in the integrand of Eq. 
is expressed as: 



cos[s\/ de 2 + be + c] + cos[sy de 2 + ee + f] 



J de 2 + ie + f] , _ „ . 
V = I ((a - d)e 2 + {b -e)e + c-f) 



cos[sy/ de 2 + be + c] — cos[sy de 2 + ee + f] 

where c = /, if we assume only colormagnetic vacuum fields. For the following, 
we make use of symbolic calculations with Mathematical of Taylor expansions 
in powers of e in particular, for which this rewriting helps to organize and cut 
down the size of the lengthy expressions. 

Thus, the ultraviolet (s — > 0) structure of the unrenormalized but QCD- 
subtracted effective Lagrangian emerges after expanding up to and including 
terms of 0(e 3 ): 

5^ r (34) 

•( tr c {[a + J][is 2 + 0(& 4 )] + 0(bes 4 ) + 0([b 2 + S 2 ] S 4 )} ) G . 

Terms which are linear or cubic in G^ v (and correspondingly in e) do not con- 
tribute here because of the Gaussian ensemble average over vacuum fields; the 
leading terms 0(e °) are cancelled by the subtraction of Cq , Eq. (p9|). Further- 
more, using Eqs. (|3l|)-([33"|), we observe that the term which does not contain 
color fields, i.e. oc s _1 er[a + d] cx s~ 1 e 2 F fa/ F ,i '' in the integrand of Eq. ([34]), 
presents the usual UV divergence which is absorbed by electromagnetic charge 
and field renormalization. Due to the ensemble average, in particular with 
(G^fjG^^c oc 5 ah {g ai gps — g a 69/3j), cf. the correlator (||), the remaining finite 
terms oc s all are proportional to e 2 Fa^F 1 ^ 1 ' in the end, besides the appropri- 
ate power of the gluon condensate (0). Therefore, they are absorbed by an 
additional finite renormalization. 
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Following the previous analysis, we obtain the renormalized QCD-modified 
Euler-Heisenberg Lagrangian: 



Ceh — —-^F^vF^" + C A — C Q 



1536Wo S 2 6 [ Ct (V5sin[ S V£]) 3 

(cos[sVI](-3[& 2 + g 2 ] + 12c[5 + d] + 26s 2 &5S + lls 2 c[6 2 + £ 2 ]) 
+ cos[3s\/S](3[6 2 + e 2 ] - 12c[5 + d] - 2s 2 6c"e + s 2 c[b 2 + e 2 ]) 
-24sVdsin[sVc\(2c[a + d\ + be)}) G , (35) 



where all charges and fields are the renormalized ones by now, in particular a = 
e 2 /47r w 1/137, and where we implemented c = / for the case of colormagnetic 
vacuum fields. Naturally, the zeroth order Lagrangian for the electromagnetic 
field appears here on the right-hand side. 

We remark that after renormalization and upon expansion of the effec- 
tive action Ceh, Eq- (|35|), the QCD effects enter at 0(e 4 ) and higher, i.e. 
affecting the nonlinear effective interaction of electromagnetic fields. In the 
following section, we calculate the first nontrivial modification of the usual 
Euler-Heisenberg Lagrangian. Some comments concerning the nonperturba- 
tive content of our results will be made shortly. 

4-2 Evaluation at 0(e A ) 

Here we expand the renormalized QCD-modified Euler-Heisenberg Lagrangian 
up to 0(a 2 ), keeping all orders in the strong coupling. We consider as an 
instructive example the case of the stochastic colormagnetic vacuum fields B a 
together with an (applied) external electric field E. Thus Eqs. ([n])-(|3^) are 
replaced by: 

a = d=^F llv F^ = -E 2 , (36) 

b=-e= -4igE ■ B a t a , (37) 
c=f = g 2 B-B . (38) 

Using these, we obtain Ceh at 0(e 4 ) from Eq. (|35|), 

EH — a r li» r 
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„ 2 f°° dx / , rr -3& 4 ab 2 a 2 . 2 . -b 4 ~ab 2 . 4 

f 2a 2 / — ( tr c c e zx { 7 H ? H - a; 2 + H x 4 

i x 3 \ ° L 1L 128c 4 165 3 8c 2j L 96c 4 245 3J 

/r -5& 4 3a6 2 a 2 . . 6 4 ab 2 h 2 . , 6 4 «, 
+ ([ 128^ + 165? " 8? ] W"T2^ + 4^ ]X " 720** } 
/r 36 4 a& 2 a 2 n 2 r -6 4 a& 2 4 . , 2 , , 



where z = m 2 /y g 2 B ■ B, and where again we rotated the contour of inte- 
gration (s — > —ix), as compared to Eq. (|35|); of course, the integral is finite 
due to the renormalization, despite the appearance of individually divergent 
contributions. 

Using the value of the gluon condensate given in Eq. (^) together with 
Eqs. (||) and the approximate up- and down-quark masses, m u « 5MeV ~ 
md/2, the parameter z of the integral in Eq. (|39| ) can be expected to fluctuate 
around a value of roughly 10 -4 and 10 -3 , respectively. The corresponding 
number for strange quarks is at least two orders of magnitude larger; their 
(and heavier quark) contributions will be completely negligible in the follow- 
ing. Furthermore, for the relevant range of small z-values, the integral is very 
strongly power dominated by the term in the integrand cx x 5 /x 3 , in the region 
where coth(x) ~ 1. Thus, we calculate: 

-i />oo / 7A \ 

4h = -4 V" + 2« 2 / dx ( tr c [720236-**] x 2 coth(x) ) 



G 



2 36 °e1£.\ ?2 

evs. \ I g 

where £ denotes the Riemann zeta function, in two arguments, and we employed 
formulae 3.5513. and 9.5211. of Rcf.Ej in the second and third equality, 
respectively. Since the color matrices involved here are diagonal, as discussed, 
the sum over the eigenvalues of the resulting matrix replaces the color trace. 

Then, recalling that the covariantly constant SU(3) colormagnetic vacuum 
fields are of the form B = ^(B 3 X 3 + B S X S ), in terms of the Gell-Mann matrices 
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^3,8, we obtain: 



, (4 ) _1 J?2 256 g a 2 / ,(E • B 8 ) 4 ^ (E • (V3B 3 ± B 8 )/2) 



C f H = ±# - ^ u \^-^ + y 



2 180V3m 6 \ |B 8 | 3 ^ \{V3B 3 ± B 8 )/2{ 




i 



(41) 



where the last equality follows from the fact that all three contributions are 
equal. This can be shown by a suitable coordinate transformation in the space 
of the colormagnetic fields, the ensemble of which is determined by the Gaus- 
sian distribution of Eq. (||) . 

Using the distribution (|J) of the stochastic fields, we finally obtain from 
Eq. (|4l]) the Q CD-modified Euler-Heisenberg effective action for the case of 
external electric fields at 0(a 2 ): 

EH J 1 2 45 v ' m e 4 

y 64V^(^:G^Gr : ) 1/2 ( ^ )4] : (42) 

45 ^— ' 5 m,- 2 rrn ' 

i—u,d, ... L 

where rrii and denote the quark masses and charges, respectively, and where 
we iacorporated the usual QED term (cx m~ ) due to e + e~ vacuum polariza- 
tion.Eiy 

For m u 5MeV ~ m^/2, q u = 2/3 = —2^, and the gluon condensate 
value of Eq. (^), the sum in the last term in Eq. ( [42] ) gives a numerical factor, 

y constant • q^m^ 6 w 3.86 (1 + 1CT 3 + ••.)•— j ( 43 ) 



where the dots indicate the negligible contribution of the heavier quarks; we 
factored out m~ A for comparison with the QED result. 

Clearly, the result of Eq. (^) is nonperturbative in the strong coupling a s 
and strongly dominated by the lightest w-quark term. 

The experimentally interesting situation with crossed electric and magnetic 
external fields, related to the vacuum birefingence mentioned in Sec. ljju can 
be studied along the same lines. Presently we presented only the simplest 
nontrivial case, in order to demonstrate that a sizeable interference between 
QCD and QED vacuum polarization appears to be possible, based on the 
present calculation. 
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A thorough investigation of this effect needs to be performed which also 
takes the space-time dependent stochastic vacuum fields intOpacpaunt, which arc 
implied by the SVM description of QCD vacuum structure.ffllij In particular, 
we recall from Sec. 1 that they are expected to fluctuate strongly on roughly the 
same length scale as the light quark loops, which we considered here. Whether 
this aspect of the nonperturbative stochastic fields tends to reduce the quark 
vacuum polarization considerably or not, still remains to be seen. 

5 Discussion 

We derived the modification of the Euler-Heisenberg effective action due to 
the nonperturbative QCD vacuum structure, which influences the vacuum po- 
larization contribution of electrically and color charged quarks. We employed 
a simplified version of the stochastic vacuum model (SVM), which has been 
successful in describing vario^JLiifrared aspects of QCD including confinement 
and bound state properties .fflEJ 

In Sec. 4.2 we obtained the light quark contribution to the effective action 
at 0(a 2 ) for the case of external (macroscopic) electric fields acting simultane- 
ously with stochastic colormagnetic vacuum fields. For simplicity, we assumed 
the QCD fiekis_to be covariantly constant, such as in earlier background field 
calculations ,E3E3 however, to be fluctuating in amplitude and (color) space di- 
rection. We found a sizeable contribution, of the same order of magnitude as 
the corresponding electron-positron term. 

The extension for the case of crossed electric and magnetic fields, such as 
employed in the experimental search for the usual QED vacuum birefringence 
effect ,0u can be obtained in the same way. 

However, we discussed and want to stress here once more the fact that in 
order to obtain quantitatively reliable results, the inhomogeneous SVM fields 
necessarily have to be incorporated. Their correlation length is expected to 
interfere in a still unknown way with the length scale of the presumably dom- 
inant light quark vacuum polarization loops. 

We conclude that a calculational scheme is necessaxy which allows to han- 
dle the important inhomogeneous field configurations au It is well known that 
applications of the Fock-Schwingcr technique are limited to very special field 
configurations, such as the covariantly constant ones considered here.Lrt3 

Furthermore, our present study indicates that also the properties of the 
quark condensate, in external electromagnetic fields in particular, may deserve 
a fresh look when nonperturbative QCD features are implemented with the 
help of the stochastic vacuum model. 

From the definition of the QCD effective action in the presence of external 
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electromagnetic fields, Eq. (Q) in Sec. 2, it follows immediately that the quark 
condensate could be calculated nonperturbatively indeed:Q 

(W) = • (44) 

dm 

Employing our Eqs. ( j30[ ) and (fi"2"l), adding up the contributions, however, the 
result is quite unsatisfactory. This must be attributed to the fact that the ho- 
mogeneous stochastic fields employed in a one-loop calculation here still present 
too crude an approximation for the relevant QCD vacuum properties. It will 
be very interesting to compare future-cesults of an improved SVM calculation 
with other phenomenological models. E3 Some related issues have already been 
discussed previously in a wider context 

Finally, we mention that the fate of nonperturbative particle production 



poles contained in the effective Lagrangian (35) of Sec. 4J |r after the stochastic 
averaging over vacuum fields also deserves further study.LrEj 
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